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The standard adiabatic approximation to phasing of gravitational waves from inspiralling com- 
pact binaries uses the post-Newtonian expansions of the binding energy and gravitational wave flux 
both truncated at the same relative post-Newtonian order. Motivated by the eventual need to go 
beyond the adiabatic approximation we must view the problem as the dynamics of the binary under 
conservative post-Newtonian forces and gravitational radiation damping. From the viewpoint of the 
dynamics of the binary, the standard approximation at leading order is equivalent to retaining the 
OPN and 2.5PN terms in the acceleration and neglecting the intervening IPN and 2PN terms. A 
complete mathematically consistent treatment of the acceleration at leading order should include all 
PN terms up to 2.5PN without any gaps. These define the standard and complete non-adiabatic ap- 
proximants respectively. We propose a new and simple complete adiabatic approximant constructed 
from the energy and flux functions. At the leading order it uses the 2PN energy function rather 
than the OPN one in the standard approximation so that in spirit it corresponds to the dynamics 
where there are no missing post-Newtonian terms in the acceleration. We compare the overlaps 
of the standard and complete adiabatic approximant templates with the exact waveform (in the 
adiabatic approximation) for a test-particle orbiting a Schwarzschild black hole. Overlaps are com- 
puted using both the white-noise spectrum and the initial LIGO noise spectrum. The complete 
adiabatic approximants lead to a remarkable improvement in the effectualness (i.e. larger overlaps 
with the exact signal) at lower PN (< 3PN) orders. However, standard adiabatic approximants of 
order > 3PN are nearly as good as the complete adiabatic approximants for the construction of 
efi^ectual templates. In general, faithfulness (i.e. smaller biases in the estimation of parameters) of 
complete approximants is also better than that of standard approximants. Standard and complete 
approximants beyond the adiabatic approximation are next studied using the Lagrangian models 
of Buonanno, Chen and Vallisneri in the test mass limit. A limited extension of the results to the 
case of comparable mass binaries is provided. In this case, standard adiabatic approximants achieve 
an effectualness of 0.965 at order 3PN. If the comparable mass case is qualitatively similar to the 
test mass case then neither the improvement of the accuracy of energy function from 3PN to 4PN 
nor the improvement of the accuracy of flux function from 3.5PN to 4PN will result in a signifi- 
cant improvement in effectualness in the comparable mass case for terrestrial laser interferometric 
gravitational wave detectors. 

PACS numbers: 04.25Nx, 04.30, 04.80.Nn, 97.60.Jd, 95.55Ym 



INTRODUCTION 



The late-time dynamics of astronomical binaries consisting of neutron stars and/or black holes is dominated 
by relativistic motion and non-linear general relativistic effects. The component bodies would be accelerated 
to velocities close to half the speed of light before they plunge towards each other, resulting in a violent event 
during which the source would be most luminous in the gravitational window. Such events are prime targets 
of interferometric gravitational wave (GW) detectors like LIGO/VIRGO/GEO/TAMA that are currently taking 
data at unprecedented sensitivity levels and bandwidths i Jj] . 

Binary coalescences are the end state of a long period of adiabatic dynamics in which the orbital frequency of 
the system changes as a result of gravitational radiation backreaction but the change in frequency per orbit is 
negligible compared to the orbital frequency itself. Indeed, the adiabatic inspiral phase is well-modelled by the 
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post-Newtonian (PN) approximation to Einstein's equations but this approximation becomes less accurate close 
to the merger phase. Additionally, there are different ways of casting the gravitational wave phasing formula - 
the formula that gives the phase of the emitted gravitational wave as a function of time and the parameters of 
the system. These different approaches make use of the post-Newtonian expansions of the binding energy and 
gravitational wave luminosity of the system^ . 



A. Standard approach to phasing formula 

The standard approach in deriving the phasing formula uses the specific gravitational binding energy E{v) (i.e. 
the binding energy per unit mass) of the system and its luminosity J^{v), both to the same relative accuracy |^. 
Including the radiation reaction at dominant order, however, is not a first order correction to the dynamics of 
the system, rather it is a correction that arises at 0[{v/c)^], where v is the post-Newtonian expansion parameter 
describing the velocity in the system and c is the speed of light ^. Thus, the phasing of the waves when translated 
to the relative motion of the bodies implies that the dynamics is described by the dominant Newtonian force and 
a correction at an order (v/c)^, but neglecting conservative force terms that occur at orders (f/c)^ and (w/c)'*. 
Such considerations have led to an approximation scheme in which one constructs the phasing of gravitational 
waves using the following ordinary, coupled differential equations: 

dip 2v^ dv T{v) 
dt TO ' dt mE'(v) ' 

where E'{v) = dE{v)/dv and to = mi + TO2 is the total mass of the binary. The phasing obtained by numerically 
solving the above set of differential equations is called the TaylorTl approximant 0. If the detector's motion can 
be neglected during the period when the wave passes through its bandwidth then the response of the interferometer 
to arbitrarily polarized waves from an inspiralling binary is given by 

hit) = — W cos[ipit) + ific], (1.2) 

where tp{t) is defined so that it is zero when the binary coalesces at time t — tc, ^pc is the phase of the signal 
at r\ = m\milw? is the symmetric mass ratio, D is the distance to the source and ^ is a numerical constant 
whose value depends on the relative orientations of the interferometer and the binary orbit. It suffices to say for 
the present purpose that for an optimally oriented source A=\. 

One can compute the Fourier transform i?(/) of the waveform given in Eq. (|1.2() using the stationary phase 
approximation: 



^ 1^ /^„-V2g.[2./tc-^C+^^(/)-V4] (1.3) 

^ D V 384 / ' ^ ' 



where the phase of the Fourier transform obeys a set of differential equations given by 

d^} dt Tim? E'(vf) 

^7 = 2^^*' 17= - o 2 J ! ■ 1-4 
df df 3vf T[vf) 

In the above expressions, including the post-Newtonian expansions of the energy and flux functions, the parameter 
Vf = {nmfy^^. The waveform Ea. (|1.3|) computed by numerically solving the differential equations Ea. H1.4|l is 
called Taylor Fl fil approximant. 

Before we proceed further, let us recall the notation used in post-Newtonian theory to identify different orders 
in the expansion. In the conservative dynamics of the binary, wherein there is no dissipation, the energy is 
expressed as a post-Newtonian expansion in e = (w/c)^, with the dominant order termed Newtonian or OPN and 
a correction at order e" — (w/c)^", n — 1,2,..., called nPN, with the dynamics involving only even powers of 
-y/e — (v/c). When dissipation is added to the dynamics, then the equation of motion will have terms of both odd 
and even powers of v/c. Thus, a correction of order (w/c)™ is termed as (to/2)PN. 



^ In the case of binaries consisting of spinning bodies in eccentric orbit one additionally requires equations describing the evolution 

of the individual spins and the orbital angular momentum, but this complication is unimportant for our purposes. 
^ Throughout this paper we use units in which G = c = 1. 
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In the adiabatic approximation of a test-particle orbiting a Schwarzschild black hole, the energyfunction E{v) 
is exactly calculable analytically, while the flux function !F{v) is exactly calculable numerically In 
addition, T{v) has been calculated analytically to 5.5PN order [TJ] by black hole perturbation theory |12| |. In 
contrast, in the case of a general binary including bodies of comparable masses^the energy function E(v) has 
been calculated recently to 3PN accuracy by a variety of metho ds 11^ IM ll^ 11^ T he flux function !F{v), 
on the other hand, has been calculated to 3.5PN accuracy pi. 120. 121I bl l23l l24 I2I IM I29IJ30I Isij up to now 
only by the multipolar-post-Minkowskian method and matching to a post- Newtonian source |32| . 

B. Complete phasing of the adiabatic inspiral: An alternative 

The gravitational wave flux arising from the lowest order quadrupole formula, that is the OPN order flux, leads to 
an acceleration of order 2.5PN\ti the equations-of-motion. This far-zone computation of the flux requires a control 
of the dynamics, or acceleration, to only Newtonian accuracy. The lowest order GW phasing in the adiabatic 
approximation uses only the leading terms in the energy (Newtonian) and flux (quadrupolar) functions. For higher 
order phasing, the energy and flux functions are retained to the same relative PN orders. For example, at 3PN 
phasing, both the energy and flux functions are given to the same relative 3PN order beyond the leading Newtonian 
order. We refer to this usual physical treatment of the phasing of GWs computed in the adiabatic approximation, 
and used in the current LIGO/VIRGO/GEO/TAMA searches for the radiation from inspiralling compact binaries, 
as the standard adiabatic approximation. We will denote the nPN standard adiabatic approximant as T{E[n],J-„), 
where [p] denotes the integer part of p. 

As a prelude to go beyond the standard adiabatic approximation, let us consider the phasing of the waves in 
terms of the equations of motion of the system. To this end, it is natural to order the PN approximation in 
terms of its dynamics or acceleration. From the viewpoint of the dynamics, the leading order standard adiabatic 
approximation is equivalent to using the OPN (corresponding to OPN conserved energy) and 2.5PN (corresponding 
to the Newtonian or OPN flux) terms in the acceleration ignoring the intervening IPN and 2PN terms. A complete, 
mathematically consistent treatment of the acceleration, however, should include all PN terms in the acceleration 
up to 2.5PN, without any gaps. We refer to the dynamics of the binary, and the resulting waveform, arising 
from the latter as the OPN complete non- adiabatic approximation. In contrast, the waveform arising from the 
former choice, with gaps in the acceleration at IPN and 2PN, is referred to as the OPN standard non-adiabatic 
approximation. Extension to higher-order phasing is obvious. At IPN the standard non-adiabatic approximation 
would involve acceleration terms at orders OPN, IPN, 2.5PN and 3.5PN, whereas the complete non-adiabatic 
approximation would additionally involve the 2PN and 3PN acceleration terms. 

Finally, we propose a simple extension of the above construction to generate a new class of approximants in 
the adiabatic regime. To understand the construction let us examine the lowest order case. Given the OPN flux 
(leading to an acceleration at 2.5PN), one can choose the energy function at 2PN (equivalent to 2PN conservative 
dynamics) instead of the standard choice OPN (equivalent to OPN or Newtonian conservative dynamics). This 
is the adiabatic analogue of the complete non-adiabatic approximant'^. Extension to higher PN orders follows 
naturally. For instance, corresponding to the flux function at IPN (1.5PN), the dissipative force is at order 3.5PN 
(4PN), and, therefore, the conservative dynamics, and the associated energy function, should be specified up to 
order 3PN (4PN). In general, given the flux at nPN-order, a corresponding complete adiabatic approximant is 
constructed by choosing the energy function at order [n-|-2.5]PN, where as mentioned before, [p] denotes the integer 
part of p. We refer to the dynamics of the binary and the resulting waveform arising from such considerations, as 
the complete adiabatic approximation. We will denote the nPN complete adiabatic approximant as T(E[n+2.5], ^n)- 

Before moving ahead the following point is worth emphasizing: The standard adiabatic phasing is, by con- 
struction, consistent in the relative PN or dcv of its constituent energy and flux functions, and thus unique in its 
ordering of the PN terms. Consequently, one can construct different inequivalent, but consistent, approximations 
as discussed in Ref. by choosing to retain the involved functions or re-expand them. The complete adiabatic 
phasing, on the other hand, is constructed so that it is consistent in spirit with the underlying dynamics, or 
acceleration, rather than with the relative PN orders of the energy and flux functions. Consequently, it has a 
unique meaning only when the associated energy and flux functions are used without any further re- expansions 
when working out the phasing formula. As a result, though the complete non-adiabatic approximant is more 
consistent than the standard non-adiabatic approximant in treating the PN accelerations, in the adiabatic case 



^ In this case one may also choose the energy function to 3PN accuracy and construct a complete approximant leading to 3PN 
acceleration 
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(a) (b) 

FIG. 1: Various T-approximants of Newton-normalized (^-derivative of) energy function E'rj-iv) / E'^{v) (left), and flux 
function Tt{v) / T'-s<i{v) (right) in the test mass limit along with the exact functions (denoted by X). Also plotted is the 
amplitude spectral density (per \/Hz) of initial LIGO noise in arbitrary units. 

there is no rigorous sense in which one can claim that either of the approximants is more consistent than the 
other. The important point, as we shall see is that, not only are the two approximants not the same but the new 
complete adiabatic approximants are closer to the exact solution than the standard adiabatic approximants. 

In our view, these new approximants should be of some interest. They are simple generalizations of the standard 
adiabatic approximants coding information of the PN dynamics beyond the standard approximation without the 
need for numerical integration of the equations of motion. They should be appropriate approximants to focus 
on when one goes beyond the adiabatic picture and investigates the differences stemming from the use of more 
complete equations of motion (see Section ^Oj. 

In the case of comparable mass binaries, the energy function is currently known up to 3PN order and hence it 
would be possible to compute the complete adiabatic phasing of the waves to only IPN order. One is thus obliged 
in practice to follow the standard adiabatic approximation to obtain the phasing up to 3.5PN order. Consequently, 
it is a relevant question to ask how 'close' are the complete and standard adiabatic approximants. The standard 
adiabatic approximation would be justified if we can verify that it produces in most cases a good lower bound to 
the mathematically consistent, but calculationally more demanding, complete adiabatic approximation. In this 
paper we compare the standard and complete models by explicitly studying their overlaps with the exact waveform 
which can be computed in the adiabatic approximation of a test mass motion in a Schwarzschild spacetime. The 
availability, in this case, of the exact (numerical) and approximate (analytical) waveforms to as high a PN order 
as {v / , allows one to investigate the issue exhaustively, and provides the main motivation for the present 
analysis. Assuming that the comparable mass case is qualitatively similar and a simple 77-distortion of the test 
mass case would then provide a plausible justification for the standard adiabatic treatment of the GW phasing 
employed in the literature'*. 

C. Non-adiabatic inspiral 

The phasing formulas derived under the various adiabatic approximation schemes assume that the orbital 
frequency changes slowly over each orbital period. In other words, the change in frequency Ai/ over one orbital 
period P is assumed to be much smaller than the orbital frequency v = P~^. Denoting by v the time-derivative 
of the frequency, the adiabatic approximation is equivalent to the assumption that /S.v — vP <^ v 01 v /v^ <^ 1. 
This assumption becomes somewhat weaker, and it is unjustified to use the approximation E = —T^ when the 
two bodies are quite close to each other. Buonanno and Damour [s^. l35l | introduced a non-adiabatic approach 
to the two-body problem called the effective one-body (EOB) approximation. In this approximation one solves 
for the relative motion of the two bodies using an effective Hamiltonian with a dissipative force put-in by hand. 
EOB allows to extend the dynamics beyond the adiabatic regime, and the last stable orbit, into the plunge phase 
of the coalescence of the two bodies [s^ IsgI IstI Is^ . 



Note, however, that th e vie w that the comparable mass case is just a r^-distortion of the test mass approximation is not universal. 
In particular, Blanchet l33l has argued that the dynamics of a binary consisting of two bodies of comparable masses is very different 
from, and possibly more accurately described by post-Newtonian expansion than, the test mass case. 
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Recently, Buonanno, Chen and Vallisneri '39| have studied a variant of the non-adiabatic model but using 
the effective Lagrangian constructed in the post-Newtonian approximation. We shall use both the standard and 
complete non-adiabatic Lagrangian models in this study and see how they converge to the exact waveform defined 
using the adiabatic approximation^. 

D. What this study is about 

In our study we will use the effectualness and faithfulness (see below) to quantify how good the various approx- 
imation schemes are. There are at least three different contexts in which one can examine the performance of an 
approximate template family relative to an exact one. Firstly, one can think of a mathematical family of approx- 
imants and examine its convergence towards some exact limit. Secondly, one can ask whether this mathematical 
family of approximants correctly represents the GWs from some physical system. Thirdly, how does this family 
of approximate templates converge to the exact solution in the sensitive bandwidth of a particular GW detector. 
In the context of GW data analysis, the third context will be relevant and studied in this paper. Although there 
is no direct application to GW data analysis, equally interesting is the mathematical question concerning the 
behavior of different approximations, and the waveforms they predict, in the strongly non- linear regime of the 
dynamics of the binary, which is also studied in this paper. The latter obviously does not require the details 
of the detector-sensitivity and it is enough to study the problem assuming a flat power spectral density (i.e. a 
white-noise background) for the detector noise. 

To summarize, our approach towards the problem will be two-pronged. First, we will study the problem as a 
general mathematical question concerning the nature of templates defined using PN approximation methods. We 
shall deal with four families of PN templates - the standard adiabatic, complete adiabatic, standard non-adiabatic 
and complete non-adiabatic (in particular, Lagrangian-based) approximants - and examine their closeness, defined 
by using effectualness and faithfulness, to the exact waveform defined in the adiabatic approximation. Since this 
issue is a general question independent of the characteristics of a particular GW detector, we first study the 
problem assuming the white-noise case. Having these results, we then proceed to see how and which of these 
results are apphcable to a specific detector, namely the initial LIGO ^. During the course of this study, we also 
attempt to assess the relative importance of improving the accuracy of the energy and flux functions by studying 
the overlaps of the PN templates constructed from different orders of energy and flux functions with the exact 
waveform. It should be kept in mind that this work is solely restricted to the inspiral part of the signal and 
neglects the plunge and quasi-normal mode ringing phases of the binary d S HI |^ ^ ,42i] . 

E. Effectualness and Faithfulness 

In order to measure the accuracy of our approximate templates we shall compute their overlap with a fiducial 
exact signal. We shall consider two types of overlaps 43, 44 4&J . The first one is the faithfulness which is the 
overlap of the approximate template with the exact signal computed by keeping the intrinsic parameters (e.g. the 
masses of the two bodies) of both the template and the signal to be the same but maximizing over the extrinsic 
(e.g. the time-of-arrival and the phase at that time) parameters. The second one is the effectualness which is 
the overlap of the approximate template with the exact signal computed by maximizing the overlap over both 
the intrinsic and extrinsic parameters. Faithfulness is a measure of how good is the template waveform in both 
detecting a signal and measuring its parameters. However, effectualness is aimed at finding whether or not an 
approximate template model is good enough in detecting a signal without reference to its use in estimating the 
parameters. As in previous studies, we take overlaps greater than 96.5% to be indicative of a good approximation. 



^ See Sec. Illll for a caveat in this approach. 

® The one-sided noise power spectral density (PSD) of the initial LIGO is given in terms of a dimensionless frequency x = f / fo 
by IS, Sh{x) = 9 X 10-'*'' [(4.49x)-5i5 + 0.16a;-* 52 + 0.52 + 0.32x2] ^^ere fo = 150 Hz and the PSD rises steeply below a lower 
cut-off fc = 40 Hz. 
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FIG. 2: Effectualness (left panels) and faithfulness (right panels) of various test mass TaylorTl and TaylorFI templates in 
detecting a signal from a (IMq, IOMq) binary in white-noise. Different lines in the panels correspond to different orders 
of the flux function. Each line shows how the overlaps are evolving as a function of the accuracy of the energy function. 
Standard adiabatic approximants T{Ey„^,Tri) are marked with thick dots. 



F. Organization of the paper 



In the next section we study the test-mass waveforms in the adiabatic approximation. We discuss the con- 
struction of the exact energy and flux functions as well as the respective T- approximants. The overlaps of various 
standard adiabatic and complete adiabatic approximants are also compared in this Section. Section lTTTl deals with 
the non-adiabatic approximation. Section ll V I explores the extension of the results in the comparable mass case. 
It presents the energy and flux functions which are the crucial inputs for the construction of the fiducial 'exact' 
waveform as well as the approximate waveforms followed by a discussion of the results. In the last section we 
summarize our main conclusions. 

One of the main conclusions of this paper is that the effectualness of the test-mass approximants significantly 
improves in the complete adiabatic approximation at PN orders below 3PN. However, standard adiabatic ap- 
proximants of order > 3PN are nearly as good as the complete adiabatic approximants for the construction of 
effectual templates. In the comparable mass case the problem can be only studied at the lowest two PN orders. 
No strong conclusions can be drawn as in the test mass case. Still, the trends indicate that the standard adiabatic 
approximation provides a good lower bound to the complete adiabatic approximation for the construction of 
both effectual and faithful templates at PN orders > 1.5PN. From the detailed study of test-mass templates we 
also conclude that, provided the comparable mass case is qualitatively similar to the test mass case, neither the 
improvement of the accuracy of energy function from 3PN to 4PN nor the improvement of the accuracy of flux 
function from 3.5PN to 4PN will result in a significant improvement in effectualness in the comparable mass case. 
As far as faithfulness is concerned, it is hard to reach any conclusion. To achieve the target sensitivity of 0.965 in 
effectualness corresponding to a 10% loss in the event-rate, standard adiabatic approximants of order 2PN and 
3PN are required for the (lOM©, \QMq) and (1.4Mo, 1.4Mo) binaries, respectively, when restricting to only the 
inspiral phase. (Be warned that this is not a good approximation in the BH-BH case since the approach to the 
plunge and merger lead to significant differences.) 



II. TEST MASS WAVEFORMS IN THE ADIABATIC APPROXIMATION 



Our objective is to compare the effectualness (i.e larger overlaps with the exact signal) and faithfulness (i.e. 
smaller bias in the estimation of parameters) of the standard adiabatic T(E^n^,J^„) and complete adiabatic 
T{E[n+2.5]T^n) approximants. As a by-product of this study, we would also like to have an understanding of the 
the relative importance of improving the accuracy of the energy function and flux function. Thus, what we will do 
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is to take all possible combinations of T-approximants^ of energy and flux functions, construct PN templates and 
calculate the overlap of these templates with the exact waveform. In all cases, the exact waveform is constructed 
by numerically integrating the phasing formula in the time-domain [TaylorTl approximant, cf. Eqs. and 
(|1.2|) ]. The waveforms (both the exact and approximate) are all terminated at v\so — I/'s/B, which corresponds 
to Fiso ~ 43 Hz for the (IMq, IOOMq) binary, Fi^^ - 86 Hz for the (1Mq,50Mq) binary and Fi^^ - 399 Hz for 
the (IMq, IQMq) binary The lower frequency cut-off of the waveforms is chosen to be Flow = 20 Hz. 

In this study, we restrict to approximants TaylorTl and TaylorFl since they do not involve any further re- 
expansion in the phasing formula and hence there is no ambiguity when we construct the phasing of the waves 
using approximants with unequal orders of the energy and flux functions. 



A. The energy function 



In the case of a test-particle m2 moving in circular orbit in the background of a Schwarzschild black hole of 
mass mi, where m2/mi 0, the energy function E{x) in terms of the invariant argument x = v"^ is given by 



1 - 2a: 
VI -3a; 



The associated w-derivative entering the phasing formula is 

dE{x) 



^exact(«) = 2w 



dx 



-rjv 



(1 _ 3^2)3/2 ■ 



(2.1) 



(2.2) 



We use the above exact E'{v) to construct the exact waveform in the test-mass case. In order to construct various 
approximate PN templates, we Taylor-expand E'^^^^^{v) and truncate it at the necessary orders. 



E'rpNiv) 



1 



81 675 19845 
" ~2 8 16 128 

137781 1852389^12 12196899 w^'' 



256 



1024 



2048 



(2.3) 



We follow 1^ in denoting the precise scheme used for constructing the approximant. 
® Here, Vi^^ is the velocity at the last stable circular orbit of Schwarzschild geometry having the same mass as the total mass mi +m2 
of the binary (we adopt units in which c = G = 1) and Fi^^ is the GW frequency corresponding to it. 
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Different T-approximants of the energy function E'rp{v) along with E'^^^^^{v) are plotted in Fig.nji. 



B. The flux function 



In the test-particle limit, the exact gravitational- wave flux has been computed numerically with good accuracy 
We will use this flux function (see Fig.^), along with the energy function given by Eq. H2.2|l . to construct an 
exact waveform in the test- mass case. In the test-particle limit, the GW flux is also known analytically to 5.5PN 
order from black hole perturbation theory [Tll | and given by 
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(2.5) 



We will use the energy and flux functions given by Eq. 1)2. 3|l - Eq. H2.5I) to construct various approximate 
templates by truncating the expansions at the necessary order. The different T-approximants of the flux function 
!Ft:{v) along with the (numerical) exact flux .Foxact(w) are plotted in Fig. 
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TABLE I: Effectualness of standard adiabatic T{E^„], Tn) and complete adiabatic r(-E[„+2.5] , ^n) templates in the test mass 
limit. Overlaps are calculated assuming a flat spectrum for the detector noise (white-noise). 
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0.9772 





8819 


0.9858 


0.8610 


0.9730 


3PN 





9971 


0.9991 


0.9677 


0.9713 





9861 


0.9946 


0.9821 


0.9886 





9965 


0.9959 


0.9756 


0.9792 


3.5PN 





9913 


0.9996 


0.9636 


0.9688 





9902 


0.9994 


0.9858 


0.9914 





9885 


1.0000 


0.9690 


0.9800 


4PN 





9937 


0.9973 


0.9643 


0.9663 





9975 


0.9996 


0.9903 


0.9914 





9968 


0.9992 


0.9769 


0.9795 


4.5PN 





9980 


0.9999 


0.9671 


0.9690 





9967 


1.0000 


0.9902 


0.9913 





9996 


1.0000 


0.9787 


0.9801 


5PN 





9968 


0.9979 


0.9661 


0.9667 





9994 


0.9994 


0.9913 


0.9914 





9992 


0.9991 


0.9790 


0.9797 



TABLE II: Faithfulness of standard adiabatic T{Eyri] , Tn) and complete adiabatic T(i?[„+2.5] j ^n) templates in the test mass 
limit. Overlaps are calculated assuming a flat spectrum for the detector noise (white-noise). 





(1A/0,1OM0) 


(lM0,5OMo) 


(1Mq,1OOA/0) 




TaylorTl 


TaylorFl 


TaylorTl 


TaylorFl 


TaylorTl 


TaylorFl 


Order (n) 




s 


C 


S 




C 




s 


C 


S 


C 




S 


C 




s 


C 


OPN 





6124 


0.1475 


0.6088 





1446 





2045 


0.4683 


0.2104 


0.4750 





2098 


0.4534 





2208 


0.4641 


IPN 





1322 


0.1433 


0.1350 





1461 





1182 


0.1446 


0.1236 


0.1508 





1395 


0.1901 





1432 


0.1994 


1.5PN 





5227 


0.4005 


0.5241 





3967 





3444 


0.3947 


0.3505 


0.3866 





3260 


0.7869 





3399 


0.7700 


2PN 





7687 


0.5707 


0.7680 





5689 





5518 


0.6871 


0.5535 


0.6827 





4377 


0.8528 





4506 


0.8486 


2.5PN 





4735 


0.5268 


0.4748 





5278 





2874 


0.3561 


0.2933 


0.3625 





2787 


0.4001 





2918 


0.4133 


3PN 





8629 


0.8165 


0.8932 





8277 





9420 


0.6317 


0.9334 


0.6222 





7579 


0.8407 





7570 


0.8194 


3.5PN 





9309 


0.9979 


0.9194 





9609 





6689 


0.9681 


0.6695 


0.9632 





5740 


0.9425 





5805 


0.9383 


4PN 





9174 


0.9303 


0.9087 





9176 





6693 


0.7227 


0.6701 


0.7230 





6129 


0.7112 





6236 


0.7159 


4.5PN 





9525 


0.9744 


0.9330 





9415 





7829 


0.9242 


0.7827 


0.9229 





7286 


0.9689 





7308 


0.9632 


5PN 





9370 


0.9392 


0.9225 





9241 





7275 


0.7417 


0.7276 


0.7420 





6972 


0.7409 





7027 


0.7500 



C. Comparison of standard and complete adiabatic approximants 

We present the results of our study in the test mass limit in four parts. In the first part we discuss our conclusions 
on the mathematical problem of the closeness of the standard adiabatic and complete adiabatic template families 
with the family of exact waveforms in the adiabatic approximation. In the next part we exhibit our results in the 
case of the initial LIGO detector. In the third part we compare the relative importance of improving the accuracy 
of the energy and flux functions. Finally, in the fourth part we compare the total number of GW cycles and the 
number of useful cycles accumulated by various standard adiabatic and complete adiabatic approximants. 

1. White-noise case 

First, we explore the general question of the closeness of the standard adiabatic and complete adiabatic tem- 
plates to the exact waveform assuming flat power spectral density for the detector noise. Figs. [5J^ show the 
effectualness and faithfulness of various PN templates for three archetypical binaries with component masses 
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Energy PN Order [n] ' 
(a) (D) 



FIG. 5: Effectualness (left panels) and faithfulness (right panels) of various test mass TaylorTl and TaylorFl templates in 
detecting a signal from a {IMq, IQMq), calculated for the initial LIGO noise PSD. Different lines in the panels correspond 
to different orders of the flux function. Each line shows how the overlaps are evolving as a function of the accuracy of the 
energy function. Standard adiabatic approximants T{Ey„^,Tri) are marked with thick dots. 

(1Mo,10Mq), (1Mq,50Mq) and (1Mo,100Mq), respectively 9. 

The central result of this study is that complete adiabatic approximants bring about a remarkable improvement 
in the effectualness for all systems at low PN orders (< 3PN). The complete adiabatic approximants converge 
to the adiabatic exact waveform at lower PN orders than the standard adiabatic approximants. This indicates 
that at these orders general relativistic corrections to the conservative dynamics of the binary are quite important 
contrary to the assumption employed in the standard post-Newtonian treatment of the phasing formula. On the 
other hand, the difference in effectualness between the standard and complete adiabatic approximants at orders 
greater than 3PN is very small. Thus, if we have a sufficiently accurate (order > 3PN) T-approximant of the 
flux function, the standard adiabatic approximation is nearly as good as the complete adiabatic approximation 
for construction of effectual templates. But at all orders the standard adiabatic approximation provides a lower 
bound to the complete adiabatic approximation for the construction of effectual templates. 

The faithfulness of both the approximants fluctuates as we go from one PN order to the next and is generally 
much smaller than our target value of 0.965. The fluctuation continues all the way up to 5PN order reflecting the 
oscillatory approach of the flux function to the exact flux function at different PN orders. It is again interesting to 
note that complete adiabatic approximants are generally more faithful than the standard adiabatic approximants. 
It is certainly worth exploring, in a future study, the anomalous cases where it performs worse than the standard. 

2. Initial LIGO noise spectrum 

Having addressed the general question concerning the closeness of standard adiabatic and complete adiabatic 
templates to the exact waveforms, we now compare the overlaps in the specific case of the initial LIGO de- 
tector. The effectualness and faithfulness of various test mass PN templates for the (1M0,1OM0) binary and 
{1Mq,50Mq) binary are plotted in Fig. I^and Fig. [HI respectively, and are shown in Tables Ulll IIVI and FVl As 
in the case of white-noise, here too we see that standard adiabatic approximants of order less than 3PN have 
considerably lower overlaps than the corresponding complete adiabatic approximants and the difference in over- 
laps between standard adiabatic and complete adiabatic approximants of order > 3PN is very small. Thus, as 
in the white-noise case, if we have a sufficiently accurate (order > 3PN) T-approximant of the flux function, the 



For the sake of convenience we also tabulate the results shown in Figs. 00 in Tables Q and ini 
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standard adiabatic approximation is nearly as good as the complete adiabatic approximation for the construction 
of effectual templates. Unlike in the white- noise case, for the initial LIGO noise spectrum the plots and Table IVl 
indicate that the faithfulness of PN templates greatly improves in a complete adiabatic treatment, for all orders 
studied. 

We also calculate the bias in the estimation of parameters while maximizing the overlaps over the intrinsic 
parameters of the binary. The (percentage) bias in the estimation of the parameter p is defined as 

(jp = X 100, (2.6) 

where Pmax is the value of the parameter p which gives the maximum overlap. Along with the maximized overlaps 
(effectualness) , the bias in the estimation of the parameters m and ry are also quoted in Tables IIIII IIVI It can be 
seen that at lower PN orders (order < 3PN) the complete adiabatic approximants show significantly lower biases. 
Even at higher PN orders complete adiabatic approximants are generally less-biased than the corresponding 
standard adiabatic approximants. 



3. Accuracy of energy function Vs. flux function 

In most of the cases, TaylorTl and TaylorFl templates show trends of smoothly increasing overlaps as the 
accuracy of the energy function is increased keeping the accuracy of the flux function constant. This is because 
the T-approximants of the energy function smoothly converge to the exact energy as we go to higher orders (see 
Fig. nji). On the other hand, if wc improve the accuracy of the flux function for a fixed order of energy, the 
overlaps do not show such a smoothly converging behavior. This can be understood in terms of the oscillatory 
nature of the T-approximants of the flux function. For example, templates constructed from IPN and 2.5PN 
flux functions can be seen to have considerably lower overlaps than the other ones. This is because of the poor 
ability of the IPN and 2.5PN T-approximants to mimic the behavior of the exact flux function (see Fig.^). This 
inadequacy of the IPN and 2.5PN T-approximants is prevalent in both test mass and comparable mass cases. 
Hence it is not a good strategy to use the T-approximants at these orders for the construction of templates. On 
the other hand, 3.5PN and 4.5PN T-approximants are greatly successful in following the exact flux function in 
the test mass case, and consequently lead to larger overlaps. 

We have found that in the test mass case if we improve the accuracy of energy function from 3PN to 4PN, 
keeping the flux function at order 3PN, the increase in effectualness (respectively, faithfulness) is ~ 0.36% (—16%). 
The same improvement in the energy function for the 3.5PN flux will produce an increase of ~ 0.36% (13%). 
On the other hand, if we improve the accuracy of flux function from 3.5PN to 4PN, keeping the energy function 
at order 3PN, the increase in effectualness (respectively, faithfulness) is ~ —0.17% (—12%). The values quoted 
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TABLE III: Effectualness of standard adiabatic T{Eyn], J~n) and complete adiabatic T{E[„^2.5], ^n) approximants in the 
test mass limit. Overlaps are calculated for the initial LIGO noise spectrum. Percentage biases a™, and in determining 
parameters m and ri are given in brackets. 





(IMclOM©) 




TaylorTl 


TaylorFl 


Order (n) 


s 


c 


s 


c 


OPN 


0.5910 (12, 5.7) 


0.9707 (36, 45) 


0.5527 (31, 28) 


0.8395 (48, 53) 


IPN 


0.5232 (22, 105) 


0.8397 (125, 69) 


0.4847 (18, 9.7) 


0.8393 (147, 74) 


1.5PN 


0.9688 (52, 51) 


0.9887 (8.3, 15) 


0.8398 (61, 57) 


0.8606 (4.7, 10) 


2PN 


0.9781 (18, 25) 


0.9942 (0.4, 0.6) 


0.8485 (32, 40) 


0.8693 (15, 22) 


2.5PN 


0.9490 (96, 68) 


0.9923 (26, 32) 


0.8963 (123, 75) 


0.9071 (49, 50) 


3PN 


0.9942 (0.3, 1.1) 


0.9989 (3.7, 6.2) 


0.8713 (16, 23) 


0.8822 (12, 18) 


3.5PN 


0.9940 (6.9, 11) 


0.9998 (0.6, 1.4) 


0.8685 (23, 31) 


0.8834 (17, 25) 


4PN 


0.9974 (6.2, 11) 


0.9996 (3.9, 6.9) 


0.8746 (23, 30) 


0.8817 (21, 28) 


4.5PN 


0.9988 (3.3, 5.5) 


1.0000 (0.8, 1.6) 


0.8795 (19, 27) 


0.8868 (18, 26) 


5PN 


0.9992 (4.0, 6.9) 


0.9997 (3.5, 5.7) 


0.8792 (21, 29) 


0.8825 (20, 28) 



TABLE IV: Same as Table IIIII except that the values corresponds to the (1M0,5OMq) binary. 



Order (n) 



OPN 

IPN 
1.5PN 

2PN 
2.5PN 

3PN 
3.5PN 

4PN 
4.5PN 

5PN 



(lMe,50MQ) 



TaylorTl 



S 



c 



0.8748 (24, 29) 
0.8101 (28, 104) 
0.9254 (21, 4.1) 
0.9610 (18, 16) 
0.9104 (21, 6.9) 
0.9968 (11, 21) 
0.9923 (13, 19) 
0.9979 (8.8, 13) 
0.9995 (7.1, 14) 
0.9994 (5.2, 7.7) 



0.9471 
0.9392 
0.9996 
0.9993 
0.9940 
0.9992 
0.9997 
0.9995 
1.0000 
0.9990 



(19, 14) 
(19, 40) 
(6.7, 20) 
(7.5, 16) 
(8.3, 0.7) 
(2.6, 10) 
(2.4, 5.2) 
(3.5, 4.3) 
(0.9, 1.9) 
(2.6, 2.4) 



TaylorFl 



C 



0.8294 
0.7662 
0.8772 
0.9113 
0.8630 
0.9500 
0.9445 
0.9560 
0.9590 
0.9634 



(21, 34) 
(23, 116) 
(18, 0.2) 
(16, 14) 
(19, 8.7) 
(11, 21) 
(12, 18) 
(8.9, 14) 
(7.0, 14) 
(5.9, 10) 



0.8974 
0.8898 
0.9590 
0.9583 
0.9574 
0.9648 
0.9679 
0.9672 
0.9698 
0.9690 



(17, 13) 
(18, 43) 
(6.4, 20) 
(7.7, 17) 
(9.1, 1.9) 
(2.7, 11) 
(2.8, 6.5) 
(3.9, 5.6) 
(1.5, 3.4) 
(3.4, 5.1) 



are calculated using the TaylorTl method for the {IMq, IOMq) binary for the initial LIGO noise PSD. The 
effectualness trends are similar in the case of the (IMq, 50Mq) binary also. // the comparable mass case is 
qualitatively similar to the test mass case, this should imply that neither the improvement in the accuracy of 
the energy function from 3PN to nor the improvement in the accuracy of the flux function from 3.5PN to 
4-PN will produce significant improvement in the effectualness in the comparable mass case. The trends in the 
faithfulness are very different for different binaries so that it is hard to make any statement about the improvement 
in faithfulness. 



4- Number of gravitational wave cycles 
The number of GW cycles accumulated by a template is defined as 0| 



1 fP^^" ]\[(F) 

Mtot = 7^(¥'iso - flow) - / dF—^ , (2.7) 
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TABLE V: Faithfulness of standard adiabatic T(i?[„] , and complete adiabatic T(-E[„+2.5] , J~n) templates in the test mass 
limit. Overlaps are calculated for the initial LIGO noise spectrum. 





(1M0,1OM0) 


(1M0,5OMq) 




TaylorTl 


TaylorFl 


TaylorTl 


TaylorFl 






s 


C 


S 




C 


S 


C 




s 


C 


OPN 





2186 


0.6272 


0.2108 





5879 


0.2134 


0.3498 





2145 


0.3593 


IPN 





1342 


0.1615 


0.1308 





1563 


0.1511 


0.2196 





1527 


0.2210 


1 ciPTvf 

1 .0-L i> 





3788 


0.4492 


0.3449 





6471 


0.2915 


0.9223 





2956 


0.9195 


2PN 





7449 


0.7633 


0.6279 





7782 


0.3613 


0.8157 





3674 


0.8318 


2.5PN 





3115 


0.3970 


0.2905 





3532 


0.2608 


0.4233 





2606 


0.4161 


3PN 





9633 


0.7566 


0.7913 





8297 


0.7194 


0.9686 





7057 


0.9323 


3.5PN 





8385 


0.9984 


0.6582 





7464 


0.4941 


0.9273 





5046 


0.9442 


4PN 





8356 


0.8909 


0.6527 





6725 


0.5960 


0.7934 





5864 


0.8131 


4.5PN 





9395 


0.9851 


0.6967 





7195 


0.7594 


0.9644 





7605 


0.9614 


5PN 





8960 


0.9129 


0.6770 





6821 


0.7344 


0.8350 





7432 


0.8579 



where ip\so and (^low are the GW phases corresponding to the last stable orbit and the low frequency cut-ofF, 
respectively, and N(F) = F"^ / F is the instantaneous number of cycles spent near some instantaneous frequency 
F (as usual, F is the time derivative of F). However, it has been noticed that ji^, the large number Atot is not 
significant because the only really useful cycles are those that contribute most to the signal-to- noise ratio (SNR). 
The number of useful cycles is defined as {44j 

A/'u.efui ^ (^f^'^° ^wif)Nif)^ ^w{f)^ , (2.8) 

where w{f) = a^{f)/h\{f). If Sn{f) is the two-sided PSD of the detector noise, hn{f) is defined by h\{f) = 
fSnif), while a{f) is defined by |-ff(/)| — a{tf)/[F{tf)]^^'^ where H{f) is the Fourier transform of the time- 
domain waveform h{t) (See Eas. H1.2|l and (|1.3|) ) and tf is the time when the instantaneous frequency F{t) reaches 
the value / of the Fourier variable. 

The total numbers of GW cycles accumulated by various standard adiabatic T{E[n],J^n) and complete adiabatic 
r(£'[„^2.5]i-^n) approximants in the test mass limit are tabulated in Table IVIl along with the number of useful 
cycles calculated for the initial LIGO noise PSD. We use Ea. H1.4() to calculate F and numerically evaluate the 
integrals in Ea. 1)2.81) to compute the number of useful cycles. In order to compute the total number of cycles, we 
numerically evaluate the integral in Ea. l2.7|l 

It can be seen that all complete adiabatic approximants accumulate fewer number of (total and useful) cycles 
than the corresponding standard adiabatic approximants. This is because the additional conservative terms in 
the complete adiabatic approximants add extra acceleration to the test mass which, in the presence of radiation 
reaction, would mean that the test body has to coalesce faster, and therefore such templates accumulate fewer 
number of cycles. Notably enough, approximants (like 3PN and 4.5PN) producing the highest overlaps with the 
exact waveform, accumulate the closest number of cycles as accumulated by the exact waveform. This is indicative 
that the phase evolution of these approximants is closer to that of the exact waveform. On the other hand, the 
fractional absolute difference in the number of cycles of the approximants producing the lowest overlaps (like 
OPN, IPN and 2.5PN) as compared to the exact waveform is the greatest, which indicates that these templates 
follow a completely different phase evolution. 

In order to illustrate the correlation between the number of (total/useful) cycles accumulated by an approximant 
and its overlap with the exact waveform, we introduce a quantity 6Afn — - — "j^^ '"-^^"^l -v^hich is the fractional 
absolute difference between the number of (total/useful) cycles accumulated by a template and the exact waveform. 
Here Afn and A/'cxact are the number of (total/useful) cycles accumulated by the nPN approximant and exact 
waveform, respectively. In Fig. \7\ we plot SAfn of various standard adiabatic and complete adiabatic approximants 
against the corresponding overlaps in the case of a (IMq, IOMq) binary. 

The following points may be noted while comparing the results quoted here for the number of cycles with those 
of other works, e.g. Refs. |25ll26Ll2^ . As emphasized in Ref. one can get very different results for the phasing 
depending on whether one consistently re-expands the constituent energy and flux functions or evaluates them 
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TABLE VI: Number of GW cycles accumulated by various standard adiabatic T{E[„j, J^n) and complete adiabatic 
T{Ey„^2.5]-,J~n) approximants in the test mass limit. The number of useful cycles calculated for the initial LIGO noise 
PSD is also shown in brackets. We choose a low frequency cut-off of 40Hz. 





(1A/0,1OM0) 


(1Mq,5OM0) 


(1A/0,1OOM0) 


Order (n) 


O 


n 


o 




n 


Q 
O 


n 


OPN 


481 (92.3) 


424 (74.6) 


118 (110) 


77.8 


(64.4) 


13.6 


6.7 


IPN 


ODU \LLt) 




1 on ( 1 Q.(<\ 

loU l^ioDJ 


1 9/1 


(iU4j 




lU.D 


1.5PN 


457 (81.7) 


433 (71.8) 


88.8 (76.3) 


58.5 


(38.2) 


8.4 


2.3 


2PN 


447 (77.7) 


440 (74.0) 


77.0 (61.8) 


62.5 


(41.5) 


6.1 


2.6 


2.5PN 


464 (84.5) 


454 (79.6) 


96.8 (85.5) 


74.5 


(50.5) 


9.7 


2.9 


3PN 


442 (74.7) 


440 (73.3) 


64.5 (45.2) 


58.1 


(35.5) 


3.4 


1.6 


3.5PN 


445 (76.1) 


442 (74.5) 


68.7 (49.7) 


60.6 


(36.8) 


4.0 


1.4 


4PN 


445 (75.8) 


443 (75.2) 


66.4 (45.1) 


62.9 


(39.0) 


2.9 


1.6 


4.5PN 


443 (75.1) 


442 (74.5) 


63.7 (42.0) 


60.0 


(35.6) 


2.5 


1.2 


5PN 


444 (75.3) 


443 (75.0) 


63.8 (40.9) 


62.2 


(37.8) 


2.1 


1.4 


Exact 


442 (74.1) 




59.1 (34.3) 






0.9 





without re-expansion. In the computation of the number of useful cycles different authors treat the function F 
differently, some re-expand and others do not, leading to differences in the results. The other important feature 
we would like to comment upon is a result that appears, at first, very counter-intuitive. It is the fact that in 
some cases the number of useful GW cycles is greater than the total number of GW cycles! A closer examination 
reveals that while for most cases of interest this does not happen, in principle its occurrence is determined by the 
ratio /r = F\o^lF\^a. To understand this in more detail let us consider the ratio M-c of the number of useful cycles 
to the total number of cycles in the case of white-noise (in a frequency band i^iow to i^iso) for which 

■ ^ - -^useful _ 5 1 - j^l^ 

Mot " 12^'- i_/y3- ^'■^> 

For /r ^ 1, Wr < 1. However, as /r increases to about /r = 0.52, Mi transits from being less than one to 
becoming greater than one! Essentially this arises due to the details of the scalings of the various quantities 
involved and the point of transition depends on the PN order and the precise form of the noise PSD. For /r ~ 1, 
the calculation of useful cycles does not make much physical sense. This explains the absence of AAusefui results 
for the (lA/0, 100A/q) binary in Table IVTl 



III. NON-ADIABATIC MODELS 



Before introducing new non-adiabatic models in this section, let us recapitulate our point of view in summary. 
Contrary to the standard adiabatic approximant which is constructed from consistent PN expansions of the energy 
and flux functions to the same relative PN order, we considered a new complete adiabatic approximant (still based 
on PN expansions of the energy and flux functions but of different PN orders) but consistent with a complete 
PN acceleration. Viewed in terms of the acceleration terms they include, the standard adiabatic approximation 
is inconsistent by neglect of some intermediate PN order terms in the acceleration. The complete adiabatic 
approximation on the other hand is constructed to consistently include all the relevant PN acceleration terms 
neglected in the associated standard approximant. These models were a prelude to phasing models constructed 
from the dynamical equations of motion considered in this section. However, we have worked solely within the 
adiabatic approximation. It is then pertinent to ask whether one can construct natural non-adiabatic extensions 
of both the standard and complete adiabatic approximants. And if so, how do their performance compare? 
Indeed, the work of Buonanno and Damour [s^ within the effective one-body approach to the dynamics did find 
differences between the adiabatic and non-adiabatic solutions. In this Section we investigate whether it is possible 
to introduce non-adiabatic formulations of the standard and complete approximants considered in the previous 
Section. 

The Lagrangian models studied by Buonanno, Chen and Vallisneri seem to be the natural candidates for 
the purpose since they are specified by the acceleration experienced by the binary system. The Lagrangian models 
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FIG. 7: The fractional absolute difference 5J\f„ between the number of cycles accumulated by various approximants and 
the exact waveform, plotted against the corresponding overlaps. Standard adiabatic T{E^„], J^n) approximants are marked 
with lighter markers and complete adiabatic T{E^„^2.5], ^n) approximants are marked with darker markers. Top panels 
show SAfn for the total number of cycles and bottom panels show SMn for the number of useful cycles. The number of 
useful cycles are calculated for the initial LIGO noise PSD and the low frequency cut-off is chosen to be 40Hz. The plotted 
results are for a {IMq, WMq) binary. 



considered in Ref. |33 can be thought of as the standard non-adiabatic approximants, since, following standard 
choices, they lead to gaps in the post-Newtonian expansion of the acceleration. Generalizing these Lagrangian 
models so that there are no missing PN terms, or gaps, in the acceleration we can construct the complete non- 
adiabatic approximants. With a non-adiabatic variant of the standard and complete approximants we can then 
compare their relative performances. However, we will be limited in this investigation because of two reasons: 
Firstly, the Lagrangian models are available only up to 3.5PN order, and higher order PN accelerations are as 
yet unavailable. Secondly, the only exact waveform we have, has however been constructed only in the adiabatic 
approximation. Even in the test mass limit, the exact waveform is not known beyond the adiabatic approximation. 
Due to lack of anything better, we continue to use the exact waveform in the adiabatic approximation to measure 
the effectualness and faithfulness of the non-adiabatic approximants. 

Thus to obtain non-adiabatic approximants, the signal is constructed by integrating the equations of motion 
directly using a Lagrangian formalism. The equations are schematically written as: 

dx dv 

For the complete non-adiabatic model of order n, all terms in the PN expansion for acceleration are retained up 
to order n + 2.5 without any gaps. For the standard non-adiabatic models, on the other hand, only terms in the 
acceleration consistent with the treatment of standard phasing are retained in the acceleration, resulting in gaps in 
the acceleration corresponding to intermediate PN terms neglected in the treatment. E.g. for n = the standard 
non-adiabatic approximation includes only the and a2.5PN while the complete non-adiabatic approximation 
would include in addition the aiPN and a2PN • Given the current status of knowledge of the two-body equations 
of motion, we have only two complete non-adiabatic approximants, at OPN and IPN retaining all acceleration 
terms up to 2.5PN and 3.5PN, respectively. The associated OPN (IPN) standard non-adiabatic approximation 
retains acceleration terms corresponding to OPN and 2.5PN (OPN, IPN, 2.5PN and 3.5PN). 
The explicit terms for accelerations for each order are given as follows [s^ \Wl li^ 

aN = 2-n, (3.2) 
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(3.6) 



(3.7) 



where n = r/r. In the above Ea. H3.5|l . the logarithmic terms present in aapN in |46| have been transformed away 
by an infinitesimal gauge transformation following [47j 

The above equations are solved numerically to attain the dynamics of the system. Then the orbital phase (j){t) 
is calculated by numerically solving the equations 



2 2 
uj r . 



(3.8) 



where the calculation of the orbital angular frequency uj assumes that the orbit is circular. Once we have the 
orbital phase, the waveform is generated using Ea. H1.2|l since the orbital phase (j) is related to the GW phase (p 
by V = 2(/). 



We thank Luc Blanchet for pointing this to us and providing us this form. 
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TABLE VII: Effectualness of the Lagrarigian templates in the test mass case for the white-noise. 





(IMq.IOMq) 


(1M0,5OM0) 


(1Mq,1OOM0) 


Order (n) 


S C 


S C 


s c 


OPN 


0.5521 0.4985 


0.5553 0.8399 


0.6775 0.7789 


IPN 


0.4415 0.4702 


0.5760 0.8327 


0.6557 0.7591 



TABLE VIII: Faithfulness of the Lagrangian templates in the test mass case for the white-noise. 





(1Mq,10Mq) 


(1M0,5OM0) 


(1M0,1OOM0) 


Order (n) 


s c 


S C 


S C 


OPN 


0.0450 0.0441 


0.1778 0.0991 


0.5959 0.1851 


IPN 


0.0471 0.0474 


0.3195 0.1235 


0.3646 0.2404 



A. Standard and complete non-adiabatic approximants in the test mass case 

In the following we discuss the results of our study for the non-adiabatic waveforms in the test mass limit. 
To determine the appropriate expressions for the acceleration in the test mass case, we start with the general 
expression for the acceleration and set 77 = in the conservative IPN, 2PN and 3PN terms. Since doing the same 
in the dissipative terms at 2.5PN and 3.5PN orders prevents the orbit from decaying, we retain terms linear in 77 
at these orders and set to zero terms of higher order in 77. In the first part we discuss results found on the problem 
of the closeness of the standard and complete non-adiabatic template families with the fiducial exact waveform. 
In the second part we extend our results to the noise spectrum expected in initial LIGO. 

1. White noise 

First, we explore the general question as to the closeness of the standard non-adiabatic and complete non- 
adiabatic templates assuming a flat power spectral density for the detector noise. Tables IVIII and IVIIII show 
the effectualness and faithfulness of Lagrangian models for the same three archetypical binaries as before: 
(IM0, IOM0), (IM0, 5OM0) and (IA/0, lOOAf0) binaries. At present, the results are available at too few PN 
orders to make statements of general trends in effectualness and faithfulness. However, the main result obtained 
for the adiabatic approximants seems to hold good again for the non-adiabatic approximants: the effectualness 
is higher for the complete non-adiabatic model as opposed to the standard non-adiabatic model. This is indica- 
tive of the fact that corrections coming from the conservative part of the dynamics (i.e. the well-known general 
relativistic effects at IPN and 2PN) make an improvement of the effectualness. However, as in the adiabatic 
case, the faithfulness of both standard and complete non-adiabatic models is very poor. But, in sharp contrast 
to the adiabatic case, here it appears that the complete non-adiabatic approximation results in a decrease in the 
faithfulness of the templates. 

2. Initial LIGO noise spectrum 

Having addressed the question concerning the closeness of standard non-adiabatic and complete non-adiabatic 
templates to exact waveforms, we now compare the overlaps in the case of the initial LIGO detector. 

Tables IIXI and |X] show the effectualness and faithfulness, respectively, of Lagrangian templates for the 
(IM0, IOM0) and (lAf0, 5OM0) binaries. In this case, we see the that the effectualness of the approximants gets 
significantly improved in the complete non-adiabatic approximation and is greater than 0.9 for all the systems 
studied in this paper. Faithfulness appears to be decreased by the use of complete non-adiabatic approximation 
(this result also is in sharp contrast with the corresponding adiabatic case where we find that the complete ap- 
proximation brings about a significant improvement in the faithfulness), but again there is no indication that 
either standard or complete templates are reliable in extracting the parameters of the system. 
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TABLE IX: Effectualness of the Lagrangian templates in the test mass case for the initial LIGO noise PSD. Percentage 



determining 


parameters m and 77 are given in 


brackets. 




{1Mq,10M@) 


(1M0,5OM0) 


Order (n) 


s c 


S C 


OPN 


0.5848 (30, 26) 0.9496 (55, 107) 


0.8741 (3.3, 9.2) 0.9835 (35, 4.2) 


IPN 


0.6762 (37, 49) 0.9273 (3.1, 27) 


0.8530 (34, 191) 0.9784 (24, 28) 



TABLE X: Faithfulness of the Lagrangian templates in the test mass case for the initial LIGO noise PSD. 





(1M0,1OM0) 


(lM0,5OAf0) 


Order (n) 


S C 


S C 


OPN 


0.2463 0.1216 


0.5048 0.1747 


IPN 


0.4393 0.1823 


0.3650 0.3119 



IV. COMPARABLE MASS WAVEFORMS 



In the case of comparable mass binaries there is no exact template available and the best we can do is to compare 
the performance of the standard adiabatic and complete adiabatic templates by studying their overlaps with 
some plausible fiducial 'exact' waveform. As in the case of the test masses, here too we will consider all possible 
combinations of the T-approximants of the energy and flux functions, construct PN templates and calculate 
the overlaps of these templates with the fiducial 'exact' waveform. In all cases, the fiducial 'exact' waveform 
is constructed by numerically integrating the phasing formula in the time-domain (TaylorTl approximant), and 
terminating the waveforms ('exact' and approximate) at viso = I/'s/B which corresponds to i^iso — 1570 Hz for 
a (1.4M0, 1.4Mq) binary and Fi^o ^ 220 Hz for a (lOM©, IOMq) binary The lower frequency cut-off of the 
waveforms is chosen to be Fiow = 40 Hz. 



A. The energy function 



Unlike in the test mass limit, the energy function E{x; rj) is not known exactly in the comparable mass case but 
only a post-Newtonian expansion, which has been computed at present up to 3PN accuracy [T^ IT3 . fislll^ IT^ . 
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reads, 
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Here also, tiiso is the velocity at the last stable circular orbit of the Schwarzschild geometry having the same mass as the total 
mass mi -|- m2 of the binary. Strictly speaking, in the comparable mass case, iiiao at nPN order should be determined by solving 
E'^{v) = where E^{v) is the ^-derivative of the nth PN order energy function. Since we found that our results are qualitatively 
independent of such considerations, we stick to the choice in the test-mass limit. 
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FIG. 8; Various T-approximants of Newton-normalized (u-derivative of) energy function E'rj^(v) / E'^{v) (left) and flux 
function / J-^^iv) (right) in the comparable mass case, along with the corresponding fiducial 'exact' functions (denoted 

by X). Also plotted is the amplitude spectral density (per VHz) of initial LIGO noise in arbitrary units. 



We use this expression truncated at the necessary orders to construct the various approximate templates. To 
compute a fiducial 'exact' waveform, we use the exact energy function in the test mass limit supplemented by 
the finite mass corrections up to 3PN in the spirit of the hybrid approximation l48l|. In other words, the fiducial 
'exact' energy E'{v\ rf) will look like 
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(4.3) 



where E'^^^^^(v) is the v-derivative of the exact energy function in the test mass limit given by Eq. (|2.2|l . The T- 
approximants of the energy function E^{v; if) as well as the fiducial exact energy E'^^^^^{v; rf) are plotted in Fig.|Sli. 
The uiso corresponding to the fiducial 'exact' energy function can be determined by solving E'^^^^^{v] ?]) = 0. This 
will yield a value ^^PN-hybnd ^ 0.4294 against the v\so — 0.4082 in the test-mass case (more precisely it is the 
^^MECO If the 77-corrections are included only up to 2PN instead of 3PN, „2PN-hybrid ^ 0.4113. It is worth 

pointing that i^fs™^^^'^" — 0.4456 a-nd it is not unreasonable to expect that, with 3PN ?7-corrections the 
differences between various different ways of determining the Iso converge. (For the purposes of our analysis, 
we have checked that there is no drastic change in our conclusions due to these differences and hence we use 
uniformly the value v\so = 0.4082). 



The flux function 



The flux function in the case of comparable masses has been calculated up to 3.5PN accuracy [l^ I20II21L |2^ 
m m El , and is given by: 



where 
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(a) 



Energy PN Order [n] 



(b) 



FIG. 9: Effectualness (Fig. a) and Faithfulness (Fig. b) of various TaylorTl and TaylorFl templates in the comparable 
mass case for the initial LIGO noise spectrum. Different lines in the panels correspond to different orders of the flux 
function. Each line shows how the overlaps are evolving as a function of the accuracy of the energy function. Standard 
adiabatic approximants T{E[„], J-'„) are marked with thick dots. Label A represents the (IOMq, IOMq) binary and label 
C represents the {IAMq, 1.4Mq) binary. 
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and the value of has been recently calculated to be —11831/9240 ~ —1.28 |23| by dimensional regularization. 

To construct our fiducial 'exact' waveform, we will use the energy function given by Eq. H4.3|) and the flux 
function 
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Y^i^kv'' +BevHnv) + {MvV + B6{r])v'^ \nv) 
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k=0 



(4.8) 



where J-cx^^ctiv) is the Newton-normalized (numerical) exact flux in the test-mass limit. The expansion coefficients 
Afc's and i?6 refer to the test- mass case and Ak{r])^s and BQ{rf) refer to the comparable mass case. The exact flux 
function is thus constructed by superposing all that we know in the test mass case from perturbation methods 
and the two body case by post-Newtonian methods. It supplements the exact flux function in the test body limit 
by all the r^-dependent corrections known up to 3.5PN order in the comparable mass case. The T-approximants 
of the flux function .Ft(i'; rf) and the flducial exact flux Tc^a-ctiv; rf) are plotted in Fig.|Sl'- 



C. Comparable mass results in the adiabatic approximation 



The effectualness and faithfulness of various PN templates in the case of comparable mass binaries are plotted in 
Fig- Eh a-nd Fig.EJ), respectively, and are tabulated in Tables IXll IXIII and IXIIII The overlaps of the fiducial 'exact' 
waveform are calculated with the TaylorTl and TaylorFl approximants using the initial LIGO noise spectrum. 
Different lines in the panels of Fig.^lcorrespond to different PN orders of the flux function. Let us note that in the 
case of comparable mass binaries the complete adiabatic approximants can be calculated, at present, at most up 
to IPN order. From the Tables IXll IXIII and lXIIII one can sec that the complete adiabatic approximation generally 
improves the effectualness of the templates at OPN and IPN orders. But, as far as faithfulness is concerned, it is 
hard to conclude that one approximation is better than the other at these PN orders. 

Even though complete adiabatic approximants are not calculated for higher PN orders, the general conclusion 
one can make from Fig.^lis that the complete adiabatic approximation of the phasing will not result in a significant 
improvement in overlaps if we have a flux function of order > 1.5PN. We, thus, conclude that, provided we have a 
sufficiently accurate (order > 1.5PN) T-approximant of the flux function, the standard adiabatic approximation 
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TABLE XI: Effectualness of standard adiabatic T{E^n], J^n) and complete adiabatic T(-E[„+2.5], -^n) approximants in the 
comparable-mass case. Overlaps are calculated for the initial LIGO noise spectrum. Percentage biases am and in 
determining parameters m and rj are given in brackets. 



Order (n) 


(1OM0,1OMq) 


TaylorTl 


TaylorFl 


S 


C 


S C 


OPN 


0.8815 (14, 0.2) 


0.9515 (3.7, 0.1) 


0.8810 (11, 0.1) 0.9486 (4.1, 0.2) 


IPN 


0.8457 (59, 0.1) 


0.8957 (45, 12) 


0.8089 (52, 0.1) 0.8671 (40, 8.9) 


1.5PN 


0.9536 (3.9, 0.3) 




0.9738 (23, 29) 


2PN 


0.9833 (0.4, 0.2) 




0.9762 (0.0, 1.9) 


2.5PN 


0.8728 (14, 0.1) 




0.8800 (9.6, 0.1) 


3PN 


0.9822 (1.5, 0.0) 




0.9738 (1.3, 0.1) 


3.5PN 


0.9843 (1.4, 0.0) 




0.9746 (1.1, 0.1) 



TABLE XIL Same as Table IxTI except that the values correspond to the (1.4Mq , I.4M0) binary. 





(1.4M0,1.4M0) 




TaylorTl 


TaylorFl 


Order (n) 


S 


C 


S C 


OPN 


0.8636 (1.4, 0.2) 


0.6993 (4.3, 0.2) 


0.8488 (1.4, 0.2) 0.7592 (4.6, 0.0) 


IPN 


0.5398 (5.0, 0.1) 


0.5639 (4.3, 0.2) 


0.5381 (5.0, 0.0) 0.5614 (4.3, 0.2) 


1.5PN 


0.9516 (0.4, 0.2) 




0.9056 (0.4, 0.2) 


2PN 


0.8751 (0.0, 0.1) 




0.9315 (0.0, 0.1) 


2.5PN 


0.8517 (0.4, 0.1) 




0.8333 (0.0, 0.1) 


3PN 


0.9955 (0.0, 0.3) 




0.9366 (0.4, 0.1) 


3.5PN 


0.9968 (0.0, 0.3) 




0.9376 (0.4, 0.1) 



provides a good lower bound to the complete adiabatic approximation for the construction of both effectual and 
faithful templates in the case of comparable mass binaries. It should be kept in mind that unlike the test mass case 
where the exact energy and flux functions are known leading to an exact waveform in the adiabatic approximation, 
in the comparable mass case we are only talking about fiducial energy and flux functions constructed from what 
is known. Probably, the fiducial waveform in this case has much less to do with the exact waveform predicted by 
general relativity. 

Tables lxJ and lxnl indicate that, to achieve the target sensitivity of 0.965 in effectualness corresponding to a 10% 
loss in the event-rate, standard adiabatic approximants of order 2PN and 3PN are required for the (lOAf0, lOM©) 
and (I.4M0, 1.4M0) binaries, respectively, when restricting to only the inspiral phase. Even though the 2PN 
standard adiabatic approximants produce the required overlaps in the case of the (IOM0, IOM0) binary, in the 
real physical case of BH-BH binaries the inspiral family would not be adequate and must be supplemented by 
the plunge part of waveform as first discussed in 6^, 35.J and later in (,35. .41^ | . A discussion of plunge requires a 
3PN description of dynamics so that the 2PN templates are no longer adequate. 

D. Comparable mass results beyond the adiabatic approximation 

Finally, for the comparable mass case, non-adiabatic waveforms were generated in the Lagrangian formalism, 
using the complete equations Eq. (|3.1I) - Eq. (|3.7|l . 

Tables IXIVI and IXVI show the effectualness and faithfulness of the standard and complete non-adiabatic La- 
grangian waveforms for the initial LIGO detector. The results are more mixed in this case than for the test 
mass case. For the OPN order, for the NS-NS binary, the standard non-adiabatic approach seems to be more 
effectual and faithful than its complete non-adiabatic counterpart. However, in the BH-BH case, the complete 
non-adiabatic seems to be more effectual; but less faithful. At the IPN order, the effectualness is always higher for 
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TABLE XIII: Faithfulness of the standard adiabatic T(_E[„],JF„) and complete adiabatic r(-E[n+2.5] , -^^i) templates in the 
comparable-mass case. The overlaps are calculated for the initial LIGO noise spectrum. 





(' 1 n A f ^ 




(1.4M0, IAMq) 




Tniilnr-T1 

1 ClyLOl 1 1 


1 ayiOl r i 


1 ay LOT 1 1 


1 ay LOT r 1 


Order (n) 


s c 


s c 


S C 


s c 


OPN 


0.5603 0.8560 


0.5608 0.8369 


0.3783 0.1624 


0.4105 0.1627 


IPN 


0.3026 0.3491 


0.3025 0.3502 


0.1520 0.1615 


0.1521 0.1614 


1.5PN 


0.7949 


0.7854 


0.7259 


0.7094 


2PN 


0.9777 


0.9601 


0.5565 


0.5777 


2.5PN 


0.5687 


0.5686 


0.5934 


0.5895 


3PN 


0.9440 


0.9446 


0.9888 


0.9246 


3.5PN 


0.9522 


0.9505 


0.9916 


0.9271 



TABLE XIV: Effectualness of the Lagrangian templates in the comparable mass case for the initial LIGO noise spectrum. 
Percentage biases Om and cr^ in determining parameters m and r) are given in brackets. 





(1.4M0, 1.4M0) 


(1OM0,1OM0) 


Order (n) 


5 C 


S C 


OPN 


0.9282 (27, 1.6) 0.5848 (32, 3.6) 


0.8533 (14, 2.5) 0.9433 (35, 8.9) 


IPN 


0.5472 (22, 1.3) 0.6439 (24, 3.1) 


0.8137 (3.6, 0.2) 0.9329 (11, 7.4) 



the complete non-adiabatic case; but faithfulness is always lower. It is interesting to note that the effectualness 
trends shown by the adiabatic and non-adiabatic approximants are the same at orders OPN and IPN. However, 
further work will be necessary to make very strong statements in this case. 

V. SUMMARY AND CONCLUSION 

The standard adiabatic approximation to the phasing of gravitational waves from inspiralling compact binaries 
is based on the post-Newtonian expansions of the binding energy and gravitational wave flux both truncated at 
the same relative post-Newtonian order. To go beyond the adiabatic approximation one must view the problem 
as the dynamics of a binary under conservative relativistic forces and gravitation radiation damping. In this 
viewpoint the standard approximation at leading order is equivalent to considering the OPN and 2.5PN terms in 
the acceleration and neglecting the intermediate IPN and 2PN terms. A complete treatment of the acceleration 
at leading order should include all PN terms up to 2.5PN. These define the standard and complete non-adiabatic 
approximants respectively. A new post-Newtonian complete adiabatic approximant based on energy and fiux 
functions is proposed. At the leading order it uses the 2PN energy function rather than the OPN one in the 
standard approximation so that heuristically, it does not miss any intermediate post-Newtonian terms in the 
acceleration. We have evaluated the performance of the standard adiabatic vis-a-vis complete adiabatic approx- 
imants, in terms of their effectualness (i.e. larger overlaps with the exact signal) and faithfulness (i.e. smaller 
bias in estimation of parameters). We restricted our study only to the inspiral part of the signal neglecting the 
plunge and quasi-normal mode ringing phases of the binary [^ l3M[3^ l40ll4]ll4 ^ . We have studied the problem 



TABLE XV: Faithfulness of the Lagrangian templates in the comparable mass case for the initial LIGO noise spectrum. 





(1.4Mq 


,1.4M0) 


(1OM0,1OMq) 


Order (n) 


S 


C 


S C 


OPN 


0.0717 


0.0658 


0.6689 0.3146 


IPN 


0.0810 


0.0771 


0.7380 0.6568 
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both for the white-noise spectrum and initial LIGO noise spectrum. 

The main resuh of this study is that the conservative corrections to the dynamics of a binary that are usually 
neglected in the standard treatment of the phasing formula are rather important at low PN orders. At the 
low PN orders, they lead to significant improvement in the overlaps between the approximate template and 
the exact waveform. In both the white-noise and initial LIGO cases we found that at low (< 3PN) PN orders 
the effectualness of the approximants significantly improves in the complete adiabatic approximation. However, 
standard adiabatic approximants of order > 3PN are nearly as good as the complete adiabatic approximants for 
the construction of effectual templates. 

In the white-noise case, the faithfulness of both the approximants fluctuates as we go from one PN order to 
the next and is generally much smaller than our target value of 0.965. The fluctuation continues all the way 
up to 5PN order probably reflecting the oscillatory approach of the flux function to the exact flux function with 
increasing PN order. Poor faithfulness also means that the parameters extracted using these approximants will 
be biased. It is again interesting to note that complete adiabatic approximants are generally more faithful than 
the standard adiabatic approximants. For the initial LIGO noise case on the other hand, the faithfulness of the 
complete adiabatic approximants is vastly better at all orders. 

To the extent possible, we also tried to investigate this problem in the case of comparable mass binaries by 
studying the overlaps of all the approximants with a fiducial 'exact' waveform using initial LIGO noise spectrum. 
It is shown that, provided we have a T-approximant of the flux function of order > 1.5PN, the standard adiabatic 
approximation provides a good lower bound to the complete adiabatic approximation for the construction of both 
effectual and faithful templates. This result is in contrast with the test mass case where we found that the complete 
adiabatic approximation brings about significant improvement in effectualness up to 2.5PN order and significant 
improvement in faithfulness at all orders. To achieve the target sensitivity of 0.965 in effectualness, standard 
adiabatic approximants of order 2PN and 3PN are required for the (lOM©, lOM©) and {IAMq, IAMq) binaries, 
respectively. Whether the complete adiabatic approximant achieves this at an earlier PN order is an interesting 
question. It is worth stressing that this result is relevant only for the family of inspiral waveforms. In the real 
physical case of BH-BH binaries the inspiral fa mily would not be a deq uate and must be supplemented by the 
plunge part of waveform as first discussed in l35| | and later in [s^, 0| ■ A discussion of plunge requires a 3PN 
description of dynamics so that the 2PN templates are no longer adequate. This is an example of the second 
variety of questions one can study in this area referred to in our introduction related to whether a template family 
indeed represents the GWs from a specific astrophysical system. 

We have also constructed both standard and complete non-adiabatic approximants using the Lagrangian models 
in Ref. |39|. However, we were limited in this investigation because of two reasons; Firstly, the Lagrangian 
models are available only up to 3.5PN order, and higher order PN accelerations are as yet unavailable which 
makes it impossible to calculate the complete non-adiabatic approximants of order > IPN. Secondly, the only 
exact waveform we have has, however been constructed only in the adiabatic approximation. So we are unable 
to make strong statements of general trends and view this effort only as a first step towards a more thorough 
investigation. From the non-adiabatic models studied, the conclusion one can draw is that while complete non- 
adiabatic approximation improves the effectualness, it results in a decrease in faithfulness. 

There is a limitation to our approach which we should point out: complete adiabatic models can be very well 
tested in the test mass where both approximate and exact expressions are available for the various quantities. 
However, complete models cannot be worked out to high orders in the comparable mass case since they need 
the energy function to be computed to 2.5PN order greater than the flux and currently the energy function is 
only known to 3PN accuracy. Also, due to the lack of an exact waveform, one is constrained to depend upon 
some fiducial exact waveform constructed from the approximants themselves. Though, in the present paper we 
have used the new approximants to construct waveform templates, one can envisage applications to discuss the 
dynamics of the binary using numerical integration of the equations of motion. 

During the course of this study, we also attempted to assess the relative importance of improving the accuracy of 
the energy function and the flux function by systematically studying the approach of the adiabatic PN templates 
constructed with different orders of the energy and the flux functions to the exact waveforms. From the study of 
test-mass templates we also conclude that, provided the comparable mass case is qualitatively similar to the test 
mass case, neither the improvement of the accuracy of energy function from 3PN to 4PN nor the improvement of 
the accuracy of flux function from 3.5PN to 4PN will result in a significant improvement in effectualness in the 
comparable mass case. 
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